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"by 
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Research  Triangle  Institute  and  University  of  Toronto 


Stmnnary 

The  distribution  of  R,  the  size  of  the  resultant  of  N  unit  vectors 
remdomly  oriented  in  three  dimensions,  was  solved  by  Rayleigh  [6]  emd  in  a  quite 
different  form  by  Fisher  [4].  These  forms  are  proved  equivalent.  An  improved 
significance  table  is  given  for  a  test  for  randomness,  proposed  by  Watson  [9], 
and  based  on  this  distribution.  An  approximate  test  is  suggested  and  compared 
with  the  exact  test.  The  power  of  the  exact  test  against  a  given  alternative 
distribution,  suggested  by  Fisher,  is  used  to  calculate  a  table  of  sample  sizes. 

1.  Introduction 


1.1  The  exact  distribution  of  R,  the  size  of  the  resultant  of  N  unit 
vectors  randomly  oriented  in  three  dimensions,  has  received  much  attention. 

It  can  be  stated  in  two  ways  which  are  not  readily  seen  to  be  equivalent. 
Rayleigh  [6 ]  gave  a  solution,  in  which  he  used  the  technique  of  Kluyver,  who 
had  first  solved  the  problem  in  two  dimensions.  Rayleigh's  density  fiinction  is 


f(R)  =  ^ 


sin  Rx  (sin  x;  dx 


(1) 


^This  work  was  done  while  the  author  was  partially  supported  by  the  Office  of 
Naval  Research. 


Page  Two 


and  he  gave  a  method  for  deriving  a  polynomial  expression  for  f(R).  This 
takes  different  forms  for  the  successive  intervals  N  >  R  >  N-2,  N-2  >  R  >  N-4, 
and  so  on  until  the  last  interval,  which  is  1>R>0  if  N  is  odd  and 
2  >  R  >  0  if  N  is  even.  Rayleigh  comments  on  the  tediousness  of  his 

procedure. 


1.2  Fisher  [k]  found  the  distribution  from  another  point  of  view.  Using  a 
geometrical  argument,  he  first  found  the  distribution  of  X,  the  component 
of  R  on  a  chosen  polar  axis.  For  three  dimensions  it  then  happens  that  the 
distribution  of  R  can  be  found  very  quickly.  Fisher's  solution  for  the 
density  function  is,  using  F(R)  to  distinguish  from  Rayleigh's  form. 


F(R)  = 


2^’^ (N-2)’- 


N/2 

^  <  N  -  R  -  2s>^"^ 


(2) 


in  which  the  notation  <x>  means  <x>  =  x  if  x  >  0 

<x>  =  0  if  X  <  0 .  (3) 

Equation  (2)  is  in  effect  a  more  compact  way  of  stating  the  polynomials 
which  arise  in  Raylei^' s  solution:  F(R)  will  be  called  Fisher's  form  of 
the  density  function.  It  forms  the  basis  of  an  exact  test  of  randomness 
discussed  in  Section  3*  i.®  easy  to  use  to  calculate  the  significance 

points  and,  with  ein  electi-onic  computer,  is  very  fast. 


1.3  A  proof  that  Rayleigh's  integral  could  be  integrated  in  the  form  of 
(2)  was  given  by  Quenouille  [5]  in  1947j  h-  method  of  induction  is  used  and 
in  fact  ths  more  general  problem  of  steps  of  different  lengths  is  solved. 

A  proof  also  appears  in  van  der  Pol  and  Bremmer  [8  ].  This  is  an  interesting 
example  of  modern  methods,  since  it  is  based  on  the  two-sided  Laplace  integral 
in  3  dimensions;  to  a  reader  unfamiliar  with  the  technique  or  the  notation, 
however,  the  proof  may  be  somewhat  difficult.  It  might  be  interesting  to 
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prove  the  equivalence  by  following  and  extending  Rayleigh's  method;  this  has 
been  done  in  Section  2.  e 


1.4  Limiting  form  as  N  becomes  large. 

The  asymptotic  distribution  as  N  becomes  large  has  many  applications. 
Rayleigh  gave  this  distribution  as 

using  the  limiting  form  of  Equation  (l).  This  distribution  is  re-derived 
from  a  statistical  point  of  view  in  Section  4,  where  it  is  used  for  an 
approximate  test  for  randomness  when  N  becomes  large.  Rayleigh's  method, 
and  other  techniques  for  random  walk  problems  will  be  found  in  a  very 
complete  paper  by  Chandrasekhar  [2  ]. 


2,  Proof  of  the  equivalence  of  Rayleigh's  and  Fisher's  forms  of  the 
distribution  function  of  R. 


2.1  It  will  be  proved -;;ic.t  --  where  f(R),  F(R)  are  Rayleigh's 

and  Fisher's  forms,  respectively  defined  in  equations  (l),  (2).  The  proof 
will  be  given  for  N  even.  It  is  required  to  prove: 


s  r  ^  ^  ax 

O 


_ 1 _ 

2^“^(N-2): 


N 


s=o 


-l)®<N-R-2s>^‘^.  (5) 


The  proof 
results: 
(a)  The 


will  follow  the  following  lines,  after  several  introductory 


(N-3) -derivative  of 
(N-3) -derivative  of 


f(R) 

R 

F(R) 

R 


is  shown  to  be  equal  to 


the 
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(■b)  Successive  integrations  are  performed,  until  the  two  forms  are  separately- 
built  up  and  shown  to  be  equal. 


2.2  Prellminaiy  results. 

It  is  of  eidvantage  to  state  two  preliminary  results, 
(a)  It  may  easily  be  shown  that,  for  N  even. 


sin  X  =  — ^  ^  -■  [2  cos  Nx 


(-1)^^^  Cy 

^s=o 


-  (^)  2  cos  (N-2)x  ...  +^|)  (-1)'^'^^  cos(N-N)x] 


(6) 


[cos(N-2s)x](-l)®  (^), 


-i(i)  • 


(b)  Expansion  of  (l-l)  gives,  for  N  even: 


n/  2  I'l/p 


s=o 


N 

N 

.2, 


(7) 


2.3  Two  Lemmas. 

Consider  the  functions 


and 


H^(R) 


H2(R) 


■/ 

o 

00 


1  ■  r,  4  N  ^ 

—  Sin  Rx  sin  x  dx 


x® 


1  T,  ^  N  ^ 

—  cos  Rx  sin  x  dx 

xP 


R)dx 


h2(x,  r)dx 


where  N  is  etn  even  integer,  p  is  an  integer  such  that  2  <  p  <  N,  and  the  variable 
R  takes  values  in  the  interval  0  <  R  <  N. 
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Lemma  1.  H^(R)  and  H2(R)  may  Le  Integrated  with  respect  to  R  by 

changing  the  order  of  integration. 

A  set  of  sufficient  conditions  for  changing  the  order  of  integration, 
and  the  formal  conclusion,  are  stated  below,  applied  to  H^(R).  (See,  e.g. , 
Titchmarsh  [7]  pp.  ^9-51 )• 

If  (a)  h^(x,  R)  is  continuous  in  the  intervals 


C  <  X  <  00 

0  <  R  <  N  , 


and 


h^(x,  r)  dR  is  uniformly  convergent  to  H^(R)  in  the 


interval  0<R<N,  then ^  dx h^(x. 


/N 

H^(R) 

o 


dR. 


These  conditions  will  now  be  examined. 

Condition  (a). 

If  0^R<N,  condition  (a)  is  easily  seen  to  hej^satisf led  for  both 

h.,(x,  R)  and  hJx,  R),  because  for  N>  p,  lim sin  Rx  and 
^  X— >0  xP 

lim  — cos  Rx  both  exist. 

X  — >0  X 

Condition  (b). 

Consider  the  integral 


sin^x 

M  (p,  N)  =  \  — ;; —  dx  for  p  an  integer  such  that  2  <  p  <  N. 


^  sir 


This  can  be  written  as 


a  .  N 


/  MS  sin^'x  ,  .  f'  sin‘’x  , 

M  (p,  N)  =  /  — r —  dx  +  /  — - —  dx. 


x^ 


N 


xP 


.  N 
sin  X 


The  first  integral  converges  since  lim 

X— >0  x^ 


exists. 
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The  second  integral  converges,  since  it  is  dominated  by  \  —  dx,  which 

converges  for  p>  1.  Thus  M(p,  N)  is  convergent.  ''l  ^  The  integrand 

of  (R)  satisfies  the  inequality 

.  .  N 

sin^x  <  for  0  <  x  <  ». 

xP  ■  xl^ 

Thus  by  the  Weierstrass  M-test,  H^(R)  is  uniformly  convergent,  since 
M(p,  R)  is  convergent.  A  similar  proof  holds  for  Hg(R).  It  follows  that 
H^(r)  and  H2(R)  are  continuous  functions  of  R. 


Lemma  2. 


H^(n)  =  I  —  (sin  x)^dx  =  0  if  p  is  odd. 


(8) 


H2(N)  =  r  ~  (sin  x)^dx  =  0  if  p  is  even.  (9) 

i  ^ 

This  may  be  shown  by  starting  with  a  result  of  Stormer  quoted  in  Bromwich 
[1],  p.  518.  This  is 


(sin  x)^dx  =  r  if  R  >  W. 

(sin\) 

If  the  integrand  is  compared  to  the  function  integral  J  may 

be  proved  uniformly  convergent  by  the  M-test.  ^  It  is  thus  continuous, 
and  satisfies  the  I’equirements  for  differentiation  under  the  integral  sign 
(Titchmarsh  [7],  p.  59)*  Hence 


-r 


sin  Ibc 
W+2. 


dJ  cos  Rx  ,  .  _  j  D  TT 

dR  "  ~  /  - iT  x)  dx  =  0  if  R  >  H. 

i 

The  equality  R  =  N  is  included  because  of  the  continuity  of  this  integral 
already  established.  Further  differentiation  gives  the  results  of  equations 
(8),  (9)  above. 
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N-3 


2. k  To  prove; 

ll\  ^“3  f(R}^  ^  F(R) 

(dRl  R  MR''.  R  ■  ■ 

Proceeding  formally,  the  left  hemd  side  of  equation  (lO)  is 


cos  Rx  /  ,  xN, 
— - -  (sin  x;  dx 


(10) 


(11) 


already  proved  uniformly  convergent,  and  the  ri^t  hand  side  is 


s=o 


Writing  (ll)  more  precisely,  and  substituting  (6)  gives 
,N-3^ 

n/2 


/d 

(dR/  \  R  / 

■  W.|  i-.  / 


dx 


s=o 
cos  Rx 


(13) 


A  typical  term  in  the  first  integral  is,  omitting  the  numerical  factor, 

/°cos  Rx  cos  (N-2s)x  , 
a 

^  -  si„2(?^)]ax  (rt) 


I  =  lim 
a—>  o 


=  lim 
a — >c 


The  integral 


lim  r 
a— >o  { 


OO  2 

'  sin  cx 


dx  =  V2j  cj  and  is  uniformly  convergent 


for  all  c . 

Hence  (l4)  becomes 


I  =  lim 


/  1  ax  -  ”/2  p!i|=S£  X  S^l 


-^o  a  X 


(15) 
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if  0  <  2s  <  N-R 
and 

poo 

I  =  11m 

a— >o  ^ 


i  dx  -  n/2  - 


N+R-2S  N-R-2S 


if  N  -  R  <  2s  <  N. 

Using  (l6),  equation  (l3)  becomes 

/d\  ^“3  /f(R)\ _ 1 

IdR;  V  R  y  -N- 


oN-2 

2  K  a 


P7n/2  \ 

- .  I(l  'i)  7 


s=o 

n/2 

s=s’  +1 


I)  *"  J" 


cos  Rx 


dx 


a 

N-R 


X 


where  s'  is  the  greatest  integer  less  than  or  equal  to  — ^ 
The  first  and  last  terms  of  (l?)  are,  using  (7) 


(-i)Vg- 


^N-1 


lim 

a  — >  o 


1  -  cos  Rx  _  ( -l) 


n/2 


X 


l) 


2 


Using  (7)  the  summation  terms  in  (17)  can  be  condensed  to 


^N-1 


n/2 


~J  R  . 


s=o 


s=o 


(16) 


(17) 


(18) 


-  (19) 
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Thus  (17)  becomes 


G4')  66  =  I  ^ 

s=o 


N> 


\R  J  2 

which  equals  the  right  hand  side  of  (12). 
Thus  using  the  <  >  notation; 


n/2 


£  r  £2^  si„“x  a-  =  ^  y  (-1)“ 

'-L  X  2  3- 


R-2s> 


(20) 


that  is, 


f(R)  F(R) 


2,5  It  is  now  possible  to  prove  the  main  resixlt: 

f(Rl-El^  . 

R  "  R 


Integrate  equation  (19)  with  respect  to  R  from  R'  to  N: 

N/_2 


sin\  dx  dR  = 


N-l 


R'  o 


2 

'^R'  s=o 


■V 


)<N-R-2s>  dR 


Applying  Lemma  1  to  the  left-hand  side  to  change  the  order  of  integration: 

.r  15/?  _  „  .  2 


(-1) 


N/2  2  1  1  cos  Rx  dR  dx 


,N 


^  ^  (.1)  (“)  <is:^ 


s=o 


When  the  R-integration  is  performed,  and  the  upper  limit  inserted,  the  left-hand 
side  vanishes  by  Lemma  2.  Thus 

m/2  _  _  _  2 


(-1) 


n/2  2 


N 


sin  R' s 


sin 


(21) 


s=o 


Dropping  the  dash  on  R,  one  may  integrate  again  with  respect  to  R  between  R'  and 
H;  and  this  process  is  continued  until  N-3  integrations  have  been  performed,  giving: 
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.  l  .w/2  2 
'  jt 


sin  Rx 
N-1 


X 


•  N  ^ 

sin  X  ox 


1 

^N-l 


s=o 


<N-R-2s> 


N-2 

(22) 


that  is: 

f(R)  _  F(R) 

R  ~  R  ■ 


Thus  the  identity  of  the  Rayleigh  and  Fisher  forms  is  established.  For 

N 

N  an  odd  integer,  the  proof  is  modified  for  the  expansion  of  sin  x  but 
follows  the  same  lines  as  that  above. 


3.  An  Exact  Test  for  Randomness. 


3.1  Since  f(R)  is  known  (Fisher's  form  will  be  used), 
the  probability  that  R  >  R^  , 


Pr  (R>Rq)  = 


f(R)dR  =a 


one  can  calculate 


(23) 


and  use  this  as  a  test  for  randomness. 

Watson  [9^  gave  this  test,  and  included  significance  points  R^  for 
a  =  5^,  and  0!  =  1^,  N  =  5  to  20.  The  values  of  R^  can  be  calculated 
by  inverse  interpolation  using  the  expression,  which  can  be  obtained  from 


(23), 


where 


Pr(R;^RQ)  =  a^(B^)  -  a^(R^)  +  a2(R^)  -  ...  , 

N-1 


(24) 


=(s)(N-i)T<  2 


si (N-s)I 


r< 


N  -  R  -  2s  '■  VNR  +N-R  -2s ^ 

o  "y  i  °  p  i 

N(N-l) 

N  _  R  _  2s  V 

.2 - >  ^NR^+N-R^-2sj  . 


(25) 


The  table  has  been  recomputed,  using  an  IBM  65O  computer  at  the  University 
of  Toronto  Computation  Centre,  and  the  values  for  four  values  of  a  are 
given  in  Table  1. 
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3.2  Exact  Test  for  Randomness. 

The  null  hypothesis  is  that  the  N  given  vectors  are  randomly 
distributed  in  3  dimensions. 

(a)  Find  R,  the  size  of  the  vector  resultant  of  the  N  vectors. 

(b)  Find  R^  from  Table  1  for  appropriate  N  and  significance 
level  0!. 

(c)  If  R>  reject  the  hypothesis  at  significance  level  0£. 
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Table  1 


w> 

4 

2.848 

3.103 

3.351 

3.490 

5 

3-186 

3.501 

3.827 

4.023 

6 

3-497 

3.853 

4.239 

4.480 

7 

3-784 

4..I78 

4.611 

4.888 

8 

4.050 

4.480 

4.955 

5.263 

9 

4.299 

4.762 

5.276 

5.612 

10 

4.535 

5.028 

5.579 

5.940 

11 

4.760 

5.281 

5.867 

6.252 

12 

4.974 

5.523 

6.140 

6.546 

13 

5-179 

5.754 

6.402 

6.832 

l4 

5.376 

5.976 

6.653 

7.103 

15 

5.567 

6.190 

6.896 

7.365 

l6 

5.751 

6.397 

7.130 

7.618 

17 

5.929 

6.598 

0  7.356 

7.863 

l8 

6.103 

6.793 

7.576 

8.100 

19 

6.271 

6.982 

7.790 

8.330 

20 

6.435 

7.166 

7.998 

8.554 

21 

6.595 

7.346 

8.200 

8.773 

22 

6.751 

7.521 

8.398 

8.986 

23 

6.904 

7.692 

8.  591 

9.194 

24 

7.053 

7.860 

8.780 

9.397 

25 

7.199 

8.023 

8.965 

9.597 

6 
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An  Approximate  Test  for  Randomness. 


4.1  Watson,  in  the  same  paper,  also  observed  that  Rayleigh's  approximation 

r2 


(Section  1.4)  for  the  density  of  R  when 
2 


W  is  large  shows  that 


N 


O 

(3.istribiited  as  X  and  hence  suggested  the  test  in  this  section.  It  might 


*130  interesting  to  derive  this  result  from  a  more  statistical  viewpoint,  which 
can  be  extended  to  the  corresponding  test  in  2  dimensions  (see,  e.g.  Curray 


(33)  and  in  fact  to  a  similar  problem  in  higher  dimensions. 

Let  X.,  y.,  z.  be  the  components  of  the  1^^  vector,  in  a  rectangular 
system  whose  origin  is  the  centre  of  the  sphere.  Using  polar  coordinates, 
z  =  cos  e,  X  =  sin  e  cos  JZ),  y  =  sin  Q  sin  j6.  The  distribution  of  has 

mean  0  and  variance 


E(z.2)  = 


2 


cos 


9  sin  0  d9  d/i)  = 


1 

3  • 


By  symmetry  the  distributions  of  y  and  z  are  the  same,  though  they  are 
not  mutually  independent. 

As  N  becomes  large,  an  application  of  the  Central  Limit  Theorem  gives 
W 

the  result  that  X  =  \  x.  will  tend  to  a  distribution  which  is  11(0,  x-)  and 

L  ^  2 

1=1 


similarly  for  Y,  Z,  and  these  distributions  will  approach  independence. 

Thus  as  the  distribution  of  are  each 

and  are  independent. 

Thus  =  I  (X^  +  +  Z^)  X  ^  approximately.  (26) 

I'J  J 

The  approximate  test  is  now  given. 


4. 2  Approximate  Test  for  Randomness. 

The  test  is  of  the  same  null  hypothesis  of  Section  3.2 

(a)  Find  R  as  in  Section  3*2. 

2  TiT  2  2 

(b)  Find  R^  from  the  equation  =  jX^  (a)  where  X^  (a)  is  the 

'^significance  ].evel  of  the  distribution,  upper  tail. 
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4. 3  Merits  of  the  Approximate  Test. 

The  true  significance  levels,  c( ,  of  values  of  R  calculated  as  in 

o 

Section  4,  have  been  calculated  using  (214-).  In  Table  2  are  given  the  correct 
values  of  R^,  the  values  of  R^  using  the  approximation  (ceilled 
the  true  level  a' ,  for  several  values  of  H. 

It  can  be  seen  that  the  approximation  gives  R^  too  large,  and  thus 
a' ^  05  as  N  becomes  larger,  however,  of — >  a. 

Table  2 


a  =  LO^ 

a  =  5^ 

a.=  2$ 

N 

R 

R  * 

cei, 

R 

R  * 

ce<^ 

R 

R  * 

Of  ^ 

0 

0 

0 

0 

0 

0 

k 

2.848 

2.887 

9.09 

3.103 

3.228 

3-^19 

3.551 

3.662 

0.4 

10 

4.535 

4.564 

9.63 

5*028 

5.104 

4.45 

5.579 

5.726 

1.5 

25 

7.200 

7.217 

9.859 

0.623 

8.070 

4.79 

8.964 

9.054 

1.82 

5.  Power  of  test  for  remdomness 

Suppose  the  test  for  randomness  is  used  when  the  vectors  in  fact  come 
from  the  distribution  suggested  by  Fisher  in  his  paper,  and  used  in 
examining  paleomagnetic  data  (see,  e.g.  Watson  and  Irving  [lOj . )  In  this 
distribution  the  density  function  of  the  spherical  polar  cooidinates  6,  ^ 
of  he  vectors  is 

g(©>  =  2  sfeh'  K. 

The  distribution  of  R  is  now 

■  *■!<.  <">  =  (nfe-R)'' 

When  fV=  0,  the  distribution  is  random;  f|^(R),  from  (27)  then  becomes 
f(R). 

The  hypothesis  of  randomness  will  now  be  rejected  with  a  probability 
X(N,  a)  given  by 
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X(N,  K,,  a) 


f|^(R)  dR 


(28) 


where  R  is  the  appropriate  significance  point  for  the  test  of  size  q;>  for 
o 

the  value  of  N,,  from  Table  1.  This  gives  the  power  of  the  test  againot  the 
above  ‘ilternatlve,  and  can  be  regarded  as  "detecting  with  probability 

X(n,  a.).  A  practical  question  is  to  ask  what  N  is  needed  to  detect  a 
given  with  a  probability  when  the  test  used  has  size  a.  I^is  has 

been  calculated  by  finding  the  value  of  .X  from  (28),  and  is  presented  in 
Table  3.  The  size  of  the  randomness  test  being  applied  is  given  first  (a), 
then  the  it  is  desired  to  detect  and  the  Table  gives  the  N  required  for 
several  values  of  X.  For  large  values  of  N,  Rayleigh's  appioximation  (^) 
has  been  used  in  (28),  both  f°i"  and  for  f|^(R)  =  this  gives  X  slightly 

too  small.  Thus  in  the  higher  values  of  the  table,  it  is  possible  that  any 
N-value  is  1  too  large,  a  conservative  error. 
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Number  N,  in  sample,  needed  to  detect  a  value  of  Yv,  with  probability 
when  applying  Watson's  test  for  randomness.  The  table  is  presented  for  k 
values  of  a,  the  size  of  the  test  of  randomness;  and  for  k  values  of  X. 
Note:  The  table  does  not  include  vaJ-ues  of  N  greater  than  70* 
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Size  of  test 

a 

Detection  prob. 

xi 

■  i 

Size  of  test 

a 

Detection  prob. 

90 

95 

98 

99 

90 

95 

98 

99 

2^ 

1.0 

56 

66 

1% 

63 

1.2 

40 

48 

56 

63 

46 

53 

62 

69 

l.h 

30 

36 

43 

48 

35 

4l 

48 

52 

1.6 

24 

29- 

36 

38 

27 

32 

38 

42 

1.8 

20 

23 

29 

31 

23 

26 

31 

3I+ 

2.0 

17 

20 

23 

26 

20 

22 

26 

30 

2.2 

15 

17 

20 

22 

17 

19 

22 

26 

13 

15 

18 

20 

15 

17 

20 

22 

2.6 

12 

14 

16 

17 

13 

15 

17 

19 

2.8 

11 

12 

14 

16 

12 

14 

16 

17 

3.0 

10 

11 

13 

14  j 

11 

13 

14 

16 

• 

3.2 

9 

10 

12 

13 

10 

12 

13 

15 

3.i^ 

9 

10 

11 

12 

10 

11 

12 

13 

3.6 

8 

9 

10 

11 

9 

10 

12 

13 

3.8 

8  ' 

9 

10 

11 

9 

10 

11 

12 

• 

7 

8 

9 

10 

8 

9 

10 

11 

4.2 

7 

8 

9 

10 

8 

9 

10 

11 

4.4 

7 

7 

8 

9 

7 

8 

9 

10 

4.6 

6 

7 

8 

9 

7 

8 

9 

10 

4.8 

6 

7 

8 

8 

7 

8 

9 

9 

5.0 

6 

7 

8 

8 

7 

7 

8 

9 
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